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Abstract 

This work is the first check of gauge invariance for nonperturbative calculations in light-front 
QED. To quantize QED in an arbitrary covariant gauge, we use a light-front analog of the equal-time 
Stueckelberg quantization. Combined with a Pauli-Villars regularization, where massive, negative- 
metric photons and fermions are included in the Lagrangian, we are then able to construct the light- 
front QED Hamiltonian and the associated mass eigenvalue problem in a Fock-space representation. 
The formalism is applied to the dressed-electron state, with a Fock-space truncation to include at 
most one photon. From this eigenstate, we compute the anomalous magnetic moment. The result 
is found to be gauge independent, to an order in a consistent with the truncation. 

PACS numbers: 12.38.Lg, 11.15.Tk, ll.lO.Gli, ll.lO.Ef 
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I. INTRODUCTION 



Any calculation in a gauge theory should be checked for its gauge dependence. Unfor- 
tunately, nonperturbative calculations in light-front QED [l| have been limited to a single 
gauge, usually light-cone gauge. This is due to the need to solve the constraint equation 
for the nondynamical part of the fermion field, which is entangled with the photon field. A 
careful use of Pauli-Villars (PV) regularization [2| has been shown jsj to allow the use of 
Feynman gauge, by providing cancellation of the photon-field dependence in the constraint 
equation. What is remarkable, however, is that this cancellation is actually not unique to 
Feynman gauge but holds for any gauge. Thus, nonperturbative calculations can be done in 
any gauge, provided the free-photon part of the Hamiltonian can be constructed. Here we 
provide such a construction for an arbitrary covariant gauge and apply the formalism to a 
calculation of the dressed-electron eigenstate and its anomalous magnetic moment, in order 
to investigate the gauge invariance of the result. 

This builds on earlier work on Yukawa theory 01 and QED fs], [sj-joj , where PV particles 
are used to regulate a light-front Hamiltonian and the eigenstates of the Hamiltonian are 
computed in one or more charge sectors of the theory. The eigenstate is expanded in a 
truncated Fock basis. The eigenvalue problem becomes a coupled set of integral equations for 
the wave functions, which are the coefficients of the Fock states in the expansion. Truncation 
keeps the coupled set finite in size, and the PV regularization keeps the integrations finite. 
For severe truncations, the coupled set can be solved analytically. In general, the set is 
solved numerically [g], [sj . The renormalization can be handled in a standard way, with the 
bare parameters of the original Lagrangian fixed by physical constraints, or by a sector- 
dependent parameterization, where the bare parameters become dependent on the Fock 
sectors connected by the terms in the Hamiltonian 0, 10-l3|- However, in a weakly coupled 



theory such as QED, such an approach cannot be expected to compete with high-order 
perturbation theory, due to numerical errors. Consideration of QED is a test for a method 
intended for strongly coupled theories. 

Light-cone coordinates [1, 15[ are used in order to have well-defined Fock-state expansions. 



which are at the heart of the method. We define these coordinates as a;"'" = t -|- ^ for time 
and X = (x~ = t — z, x±) for space, with x± = {x, y). The light-cone energy is p" = E — pz 
and momentum, p = (p+ = E +Pz,P±)- The mass-shell condition p"^ = m? relates these as 
p~ = {rn? +p\)/p'^. The positivity of p+ keeps the vacuum simple and prevents vacuum 



contributions to the Fock expansions, except for the possibility of zero modes p^. These 
modes of zero p"*" can be neglected in theories where symmetry breaking does not occur. 

Our new construction of the Hamiltonian is based on a light-front analog of the equal-time 
Stueckelberg quantization of a massive vector field [l^S The Stueckelberg quantization is 
known to allow for a zero- mass limit. It is also useful as a way to treat the physical and PV 
photons on an equal footing, consistent with the need to maintain the PV regularization. The 
quantization adds a fourth (unphysical) polarization to the three physical polarizations. The 
unphysical polarization is the only one that does not satisfy the Lorentz gauge condition 
d ■ A = However, it does satisfy the Euler-Lagrange field equation, because its four- 
momentum is placed on a different, gauge-dependent mass shell, chosen in just such a way 
as to satisfy the field equation. The key to the light-front analog is that the chosen mass 
shell is invoked for the minus component of the momentum, rather than the zero component. 

^ For an alternative construction for the massless case, which uses the canonical Dirac constraint procedure, 
see [18|. 
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The details of this can be found below, in Sec. HIl 

With this new quantization, we can formulate mass eigenvalue problems for the eigen- 
states of QED in an arbitrary covariant gauge and test for the gauge invariance of physical 
quantities computed from the eigenstates. We expect that gauge invariance will be bro- 
ken by approximations made in solving the eigenproblems. One such approximation is the 
Fock-space truncation used to reduce the eigenproblem to a finite size. Another is retention 
of finite values for the regulating PV masses; the regularization is constructed with use of 
flavor- changing currents that explicitly break gauge invariance [3| and are removed only in 
the infinite-PV-mass limit, which may not be possible in a numerical calculation. 

As a first test, we apply this formalism to a calculation of the dressed-electron eigenstate. 
Fock space is truncated to include only the bare-electron state and the one-electron/one- 
photon states, plus their PV analogs. This leads to an analytically solvable problem, reduced 
to an effective 2x2 matrix problem in the one-electron sector. From the solution, the 
anomalous moment can be computed, from the zero-momentum-transfer limit of the spin- 
flip transition amplitude. 

To obtain meaningful results, it is important to maintain the chiral symmetry of the 
massless-electron limit. This is achieved by adjusting the coupling strengths of the PV 
photons, to ensure that the dressed mass is zero when the bare mass is zero. As will be 
seen below, in Sec. IIIIl this requires two PV-photon flavors in any gauge. However, one 
flavor is sufficient if the limit of infinite PV-fermion mass is taken. This also holds for 



the sector-dependent approach 13|], as we show in Sec. IIIIl In general, the constraint of 
chiral-symmetry restoration and the PV-photon couplings are gauge dependent. 

There are, of course, other nonperturbative methods. Lattice gauge theory [l9| is partic- 
ularly successful, and use of Dyson-Schwinger equations [2(1 has produced notable results. 
However, these lack the direct access to wave functions in Minkowski space, which light-front 
Hamiltonian methods provide [l| . Thus, the methods are quite complementary, particularly 
now that light-front calculations can be done in an arbitrary gauge. There are also a light- 



front lattice method, the transverse lattice j2l|; a light-front approach in terms of effective 



fields [22J; and a supersymmetric formulation for discrete light-front Hamiltonians specifi- 



cally for supersymmetric theories [23 



The remainder of the paper contains the following sections. The general formalism for an 
arbitrary covariant gauge is given in Sec. [Tll The dressed-electron eigenproblem is solved in 
Sec, mil and used there to compute the anomalous moment. Section HVl provides a summary 
of the method and of the results obtained. Some details are left to three Appendices. 



II. LIGHT-FRONT QED IN AN ARBITRARY COVARIANT GAUGE 

We begin with the QED Lagrangian for Lorentz gauge with an arbitrary gauge parameter 
( and additional PV fields: 

2 

2 
i=0 



(2.1) 



3 



Here 

2 2 

^ = Zlv^^*' = V^^^A^' Fif^u = d^Ai^ - d^A^,, (2.2) 

i=0 i=0 

i = corresponds to a physical field, and i = 1 and 2, to PV fields. The photon fields have 
mass /ij, and the zero-mass limit /io — )■ for the physical photon is to be taken later. 

The coupling coefficients /3j and satisfy constraints. To keep e as the charge of physical 
fermion, we set /3o = 1 and = 1- To regulate ultraviolet divergences that come from 
loop integrals, we arrange cancellations for each internal line summed over physical and PV 
fields, by imposing the constraints 



0. 



(2.3) 



i=0 



i=0 



Two remaining coefficients, say ^2 and /32, are fixed by requiring chiral-symmetry restoration 
in the massless-electron limit jsj and a zero photon eigenmass [9|. 

The dynamical fields are and Ai^. We quantize the dynamical fermion fields in the 
usual way 

1 



V16 



with 



{hsik),bl,^,ik:} = (-lysu'Sss'Sik-k!) 



For the vector fields, we apply a light-front analog of Stueckelberg quantization [17 
Consider the Lagrangian of a free massive vector field: 



-\F' + \l^A--\ad-Af. 



The Euler-Lagrange field equation is 



[u + ^^)A,-{l-Od,{d-A) = Q. 



This equation is satisfied by the Fourier expansion 

3 



A^{x) 



dk 



. A=l 



+e 



^\k) 



ao{k)e-'^'' + al{k)e'^'' 

with k a four- vector associated with a different mass /t = fi/^/C, such that 

~k = k, k- = {kl + fi^)/k+. 
The polarization vectors are defined by 

e(i'2)(fc) = (o,2ei,2-fc±A+,ei,2), 



;(3) 

,(0) 



ik) = iikl-^')/k^,k^,k^)/fx, 



e^"Kk) = k/ii, 



(2.4) 



(2.5) 
(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 

(2.11) 
(2.12) 
(2.13) 
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and satisfy k ■ e^^^ = and e'^'^^ ■ e*^'*''^ = —Sw for A, A' = 1, 2, 3. The first term in satisfies 
(□ + /i^)^^ = and d ■ A = separately. The A = term violates each, but the field 
equation is satisfied. The gauge condition 5 ■ A = is to be satisfied by projection of states 
onto a physical subspace. 

The light-front Hamiltonian density is 



n = H|c=i + -{1-C){d- A){d- A- 2d^A- - 2d± ■ Aj 



(2.14) 



with the Feynman-gauge piece being 



^J.=o 



The metric of each field component is defined by e'^ 
tonian for the free massive field is then found to be 



(2.15) 



-1, 1, 1, 1). The light-front Hamil- 



dxH\x^ 



k+ 



(2.16) 



with = fj, for X = 1, 2, 3, but /io = yU = /i/ VC- The nonzero commutators are 

[ax{k),a[,{k:)] = e%y5{k-k:). 



(2.17) 



Thus, the Hamiltonian for the free photon field takes the usual form except that the mass 
of the fourth polarization is different and gauge dependent and that the metric of this 
polarization is opposite that of the other polarizations. In Feynman gauge, this reduces to 
the usual Gupta-Bleuler quantization 2^ . 

The nondynamical components of the fermion fields satisfy the constraints (i = 0, 1, 2) 



(2.18) 



(^tV) 



Ayma°'tpi+. 



Ordinarily, light-cone gauge {A_ = 0) is chosen, to make the constraint explicitly invert- 
ible. However, the interaction Lagrangian has been arranged in just such a way that the 
A-dependent terms can be canceled between the three constraints j^3l]. Multiplication by 
(— 1)*a/^ and a sum over i yields 



(2.19) 



as the constraint for the composite field that appears in the interaction Lagrangian. This 
constraint is the same as the free- fermion constraint, in any gauge, and the interaction 
Hamiltonian can be constructed from the free-field solution. 

Without this cancellation of A-dependent terms, the constraint would generate four-point 
interactions between fermion and photon fields, the instantaneous- fermion interactions [l[. 
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FIG. 1. The infinite-PV-mass limit of a tree graph involving an intermediate PV fermion of mass 
mi. The external fermions are the physical ones, with mass niQ. The limit yields an instantaneous 
interaction, denoted by the bar through the intermediate line. 



The addition of the PV- fermion fields has, in effect, factorized these interactions into flavor- 
changing photon emission and absorption three-point vertices. The instantaneous interac- 
tions are recovered in the limit of infinite PV fermion masses, because the light-cone energy 
denominator with an intermediate PV fermion cancels the PV-mass factors in the emission 
and absorption vertices, as illustrated in Fig. [T] 

The light-front Hamiltonian, without antifermion terms, is 



V- = J2 dp-^^{-iybUp)ks{p) (2.20) 

is ^ P 

IX 
ijlsX 

+ 6,,_.(g)f/^. g)] 4^\q - Mm -p) + H-c] • 



The instantaneous-photon terms associated with light-cone gauge do not appear. The po- 
larization vectors e^'*''' have an additional flavor index /, because they depend on the mass of 
the photon flavor. The vertex functions are given by [3| 



VUbq) = 


e p± 


■ q± ± ip± X + miTTij + p~^q~^ 




p^q^Vq'^ — J9+ 


v^j±ip,q) = 


-e p± 


■ q± ± ip± X + mirrij — p^q^ 




p^q^Vq~^ — p+ 


Vi^{p,q) = 


e p+ 


(q^ ± iq"^) + q^(p^ =1= ip'^) 




p~^q~^Vq^ — p+ 


KUp.q) = 


e 


(g^ Tiq^) + q'^ip'^ ± ip^) 




P'^q'^Vq'^ — p+ 



(2.21) 
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and 



±e mj{p^ ± ip'^) — mi{q^ ± zg 
Vl67r3 p+q+y/q+'^JF 
±e mjg+ — nijp^ 

V 167r3 p+q+^/q+ - p+ ' 

+ 



2^ 



Uij±{p,q) = r-— , , , , ' , • 

The extension to include antifermion terms is straightforward joj. We now apply this for- 
malism to a nonperturbative calculation of the dressed-electron state and its anomalous 
magnetic moment. 



III. DRESSED-ELECTRON EIGENSTATE 
A. Eigenvalue problem 

We wish to solve the light-front eigenvalue problem V~\i\){Py) = ^l^iE))- The eigen- 
state \ip{P_)) is expanded in a Fock basis where is diagonal and P± is zero. Here 
we consider the lowest order truncation, to include the bare-electron state and the one- 
electron/one-photon states. For a total = the eigenstate is of the form 

IV'^(P)) = J2ztbl{P)\0) + Y.j dhCfts^k)hUP-k)a],{m). (3.1) 

i ijsX 

The normalization condition is 

{r'{E)\r{p)) = KE-m'^- (3.2) 

For this truncation, we can remove the second PV-fermion flavor, since it plays no role in 
the regularization or chiral-symmetry restoration [5;]. We let 1712 — )■ 00, /32 — > 0, and /3i — )■ 1. 

The eigenvalue problem for this state reduces to a set of coupled equations for the bare- 
electron amplitudes and the two-body wave functions: 



X 



[v^:^{p - k, p)cf^{k) + u^:^{p - k, p)c},%{k)] 



and 



,r- - !!i±^i - di±m c}-ik) (3.4) 

l-y y \ 

IiJ2^-iy'z^P^e^\k) [V^^,^{P~k,m± + U;'±iP~k,P)Ss^] , 



7 



with y = /P~^ being the photon's longitudinal momentum fraction. The second equation 
(13. 4p is trivially inverted to find the two-body wave functions: 

(3.5) 



i' 



J^2 _ "^'+^1 _ M?A+fcj 

1-3/ y 

Substitution into the first equation (13. 3p yields a 2x2 matrix eigenvalue problem for the 
one-body amplitudes zf: 

(M^ - m2)z± = 2e^ ^(-l)^'^ [J + A J + m,mi,(Jo + A/q) (3.6) 



with 



JiM"^) 



-2(m, + m,0(/i + A/i)] 

dydkl ^ {-ly+'ii 



167r2 ^ ^2 _ !!!i±^ _ /f?±^y(l -1/)" 

J' l-y y 

+ k]_ 



dydk"] 



167r2 



E 



i-y y 



M2 - - t^±± y(l - 2/) 

i-y y 



2' 



and the gauge-dependent parts 



dydk'j 



(3.7) 
(3.8) 

(3.9) 



X 



m^y^ + k^ 



M2 - 



M2 - 



1-?/ 



j+imj^i f dydk 



1/2(1 -?/)2 



(3.10) 



X 



M2?/(1 - y) - rri^.y - k^ 



M2 



M2 ■ 
dydk'j 



(3.11) 



X 



-^^) (l-y)^ + m2fc2 



^2 _ "^j+fej _ 



2/ 



^2 _ "'+^1 _ 

l-'i/ 1/ 



For AJ we have taken advantage of the fact that ^ji{~'^y^ c,i = to simplify the expression 
by elimination from the numerator terms that are proportional to the denominator. 

For /io < M — mo, there is a line of poles in the {y, k\) plane, in an arc between the points 

aty = y±= {M^ - + /ig) ± ^{M^ - mg + /ig)2 - AM'^/j.^ /{2M^) on the longitudinal 

axis. For hq = 0, these reduce to ?/_ = and y+ = 1 — ml/M'^, as considered previously jsj. 
As in that case, we define integrals with these poles as principal values. Also, the same 
considerations hold for poles associated with denominators containing /io. 
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B. Analytic solution 

The matrix problem (13. 6p can be solved analytically, in terms of the defined integrals. 
The solution is facilitated by the identity J + A J = M'^[Iq + A/q), which was shown for 
Feynman gauge in [5| and is extended to an arbitrary gauge in Appendix El 

The analytic solutions are 

_ (M±mo)(M±mi) 

" 87r(mi - mo) [2{h + A/i) ± M(/o + A/q)] ' ^^'^^^ 

with 

^1 = TTT ^0, 3.13 

M ±mi 

and M = rrie, the physical electron mass. The solution with the lower sign is the physical 
one, because M = mo when a_ = 0. 

We fix ,^2 by requiring chiral-symmetry restoration, that is, M = for mo = when q;_ 
is equal to the physical value of a. This implies that (Ji + A/i) must be zero. From earlier 
work [S], we know that 

imu^o = i^E(-i)'^'r^^i-(/^?/-i)' (3-14) 

so that we only need to evaluate 

"^i(l-C)v^/ f {mly + k]_)dydk\ 



+ -y) + kl){mly + jj^lH - y) + kl) ' 

(3.15) 



Next, we write k^ in the numerator as k^ + m\ — m\ and use the technique in Appendix lAl 
to conclude that the fc^ + mf combination integrates to zeroJl After separation of the 
denominator into two terms, we can easily perform the remaining k\ integration, to obtain 

All Uo=o = 77-^ — / dy\n 2 , -2/1 T • 3.16 

The integral over y yields 

A/,(o)U„.o = ^ l^(-i) cK-/.f)K-/.f/C) • ^'-''^ 

Therefore, the constraint from chiral-symmetry restoration is 

(3.18) 

Thus, in any covariant gauge, two PV-photon fiavors are required to maintain the chiral 
limit. 



^ In general, k\ + mf is replaced by M^(l — y)^, but here M — 0. 
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In the limit of infinite mass for the PV electron, and with use of ^;(— = O5 the 
general constraint reduces to 



C-1 

c 



^(-l)'61n(^f//^^) = 0, 



(3.19) 



with /i any mass scale. Given ^0 = 1 and = 1 + C,2, this is solved by 

ln(/io//ii) 



6 



ln(/i2//ii 



> 0. 



(3.20) 



If the limit ^2 ^ 00 can be taken, ^2 is reduced to zero and the second PV photon flavor is 
removed. 

For the sector-dependent approach, the analogous quantity to consider in Yukawa theory 
is A{M'^), defined in Eq. (Bl) of [l^]- Combination of various pieces used in this definition 
yields, in the notation used there. 



1 f dxdR\ sr-^ 



M2 



(3.21) 

Comparison with our ([32D shows that A{M'^) = \li{M^)\^2=o- When M = and mn = 0, 
the i = 1 term of A is not zero, unless mi — ?■ 00. Therefore, Sm2 in (29b) of [13| is 
also not zero, in contradiction of chiral-symmetry restoration. Thus, the sector-dependent 
approach requires a second PV-photon flavor in Yukawa theory. For QED, the situation is 
not materially different. 

To complete the analysis of the eigensolution, we need to consider the gauge dependence 
and infrared dependence of the integral combination 



A = 2A/i(M2) - MAIo{M^ 



(3.22) 



which enters the denominator of fl3.12p . The combination 2/i(M^) — MJo(M^) also appears 
there, but is gauge-independent by definition and is known to be infrared safe jsj. From the 
definitions ( 13. 9 P of AJq and A/i, we can obtain, after eliminating k"^ from the numerator in 
the same manner as before, 



A 



1-C 

327r2C 



dydkj 



(M - mj)2[M(l -y) + rrij 



The integral yields 



A 



327r 



1 v^. 



_1)J+'- (M-m,-)^ 



(3.23) 
(3.24) 



X / dy[M{l-y) + m-\\n 



m% + fL^{l-y)-M^y{l-y) 



mh + /if (1 -y)- M^y{l - y) 



The j = term is of order (M — mof' oc a^, and the j = 1 term is of order 1/mi, after 
invocation of the chiral constraint fl3.18p to eliminate the leading, order-mi term. Thus, A 
breaks gauge invariance only in ways to be expected; gauge invariance can be attained only 
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without truncations and then only in the mi — )■ oo hmit. The high-order a correction is a 
signal of a truncation effect, and, of course, the 1/mi contribution disappears as mi — )■ oo. 

To study the dependence on the IR mass scale /iq, we consider the /xq — ^ limit of the 
j = I = term in f l3.24p . which is 

The behavior near y = is such that the term has a log divergence multiplied by {( — 1){M — 
mo). Thus, this contribution is of order a, not a^, and, of course, is absent in Feynman 
gauge. The order of the contribution is, however, still consistent with being a truncation 
error. Also, the presence of an IR divergence is consistent with the presence of other, UV 
divergences, which are uncanceled due to truncation, as discussed in j3] and as discussed 
below, with respect to normalization of the eigenstate. 



C. Anomalous magnetic moment 

We compute the anomalous magnetic moment of the dressed electron from the spin-flip 
matrix element of the electromagnetic current J j25|. The plus component of the current 
is used because, in the absence of vacuum polarization, it is not renormalized js], |2^. In 
general, the transition amplitude for absorption of a photon of momentum g by a dressed 
electron is given by 

{r{P + q)\^\^HP)) = 25.±Fi(g^) ± ^-^^6.^F,{q% (3.26) 
where Fi and F2 are the usual Dirac and Pauli form factors. The anomalous moment is 



Qe = i^2(0); normalization of the state is equivalent to Fi{0) = 1. As described in [25|, the 
limit of zero momentum transfer for F2 can be written as 

ae = TMj2f dke\-ir^C^t'{k)y ± CfAk)- (3-27) 

jls\ ' 

This form assumes complete separation of the internal and external momentum variables in 
the wave functions C*j7^, which does occur for components in terms of polarizations. The 
sum over polarizations A does not include the gauge projection, because gauge invariance 
has already been broken by both the truncation and the flavor-changing currents. The 
normalization condition (13. 2p . or -Fi(O) = 1, becomes 

1 = (4)' - + E / dke\-\r'\Cf^m\ (3.28) 

jls\ 

which determines . 

The reduction of the expression for the anomalous moment is given in Appendix [Bl The 
result in the limit of infinite PV electron mass is given in (IB7p . The normalization condition 
is evaluated in Appendix [Cl with the same infinite mass limit yielding f ICSP . From these 
expressions, the anomalous moment can be computed for various values of the UV scale /ii, 
the IR scale and the gauge parameter C^. Sample results are given in Figs. [21 [3l and 
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FIG. 2. The anomalous magnetic moment Og for the dressed-electron state truncated to include at 
most the one-electron/one-photon Fock states, as a function of the PV-photon mass /ii. The IR 
mass scale fiQ is 0.001 me, and the gauge parameter C is 1/2, 1, and 10. 
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FIG. 3. Same as Fig. [2] but as a function of the IR mass scale /hQ) with ni = 200 me 



m Because truncation errors do not allow either fiQ — )■ 0, except in Feynman gauge, or 
/ii — 7- oo, we look for regions in each where the physical quantity is relatively fiat. In Fig. HJ 
we then look for sensitivity to the gauge parameter when /iq and fii have values in such 
regions. The plot shows little sensitivity and, therefore, approximate gauge independence, 
except for small values of (. For small (, the theory is near the singular limit where the 
gauge-fixing term is removed from the Lagrangian (12. ip and truncation errors are amplified. 
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FIG. 4. Same as Fig. [2] but as a function of the gauge parameter with fiQ = 0.001 m-e and 
m = 200 me- 

The remaining gauge dependence can be seen to be consistent with the order of the 
truncation in the calculation. Because M — tjiq is of order a, we have, for the leading term 
in Oe, as given by flB7p . 



Up to normalization, the Schwinger result [27'] of a;/27r is recovered, and the gauge-dependent 
contributions, as well as some physical nonperturbative contributions, are higher order in a, 
consistent with the truncation to one photon in the Fock basis. 

To complete the analysis, we must consider the normalization factor zq, as given in (ICSp . 
At the present order of truncation, a value other than 1 for zq represents a truncation error, 
in the sense that contributions occurring at the same order in a for the numerator of the 
expectation value for the anomalous moment have been left out by the truncation; therefore, 
the gauge dependence of Zq must be due to truncation errors. 

In the IR limit /iq — )■ 0, the normalization zq does have a singular contribution of the 
form 3^i^/^. This comes from a combination of the last two terms of the curly bracket 
in (1C8|) for /iq — )■ and / = 0. Therefore, the normalization contains an IR divergence in 
addition to its usual UV divergence, except in Feynman gauge, where there is only a UV 
divergence. 

These divergences are the characteristic "uncanceled divergences" caused by Fock-space 
truncation [3|. They arise in both the standard and sector-dependent parameterizations, 
although in the latter case the IR divergence is present even in Feynman gauge [7|. For the 
standard parameterization used here, we find the divergences in the normalization factor. 




(3.29) 



and, therefore. 



«e = — 2o + C(a^l^?)• 



(3.30) 
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For the sector-dependent parameterization, the divergence is most easily seen by considering 
the probabihty for the one-electron/one-photon sector. This should be between zero and 
one, but the renormalization of the sector-dependent coupling absorbs the divergence in the 
normalization factor 1/zq and allows the probability of the one-electron/one-photon sector 
to diverge. In this case, the overall norm of the eigenstate is maintained only because the 
probability for the bare-electron sector goes to negative infinity in such a way that the sum 
of probabilities is formally one. 

Because of the uncanceled divergences, not all of the PV masses can be taken to infinity 
and, except for our standard parameterization in Feynman gauge, the physical photon mass 
cannot be taken to zero. As argued in jsj, the errors introduced by these limitations are to 
be minimized by seeking ranges of mass values over which results do not change significantly. 
This strikes a balance between the errors caused by the presence of unphysical PV fields and 
a nonzero photon mass and the errors associated with Fock-space truncation. The former 
decrease with increasing PV masses (as the PV fields are removed from the spectrum) and 
decreasing photon mass; the latter, the truncation errors, increase with increasing PV masses 
as the uncanceled divergences assert themselves. 



IV. SUMMARY 



We have developed a formalism whereby nonperturbative calculations can be done for 
light-front QED in an arbitrary covariant gauge. The formalism combines a light-front 
Stueckelberg quantization for the free photon field with a Pauli-Villars regularization that 
simplifies the constraint equation for the nondynamical part of the fermion field. The Stueck- 
elberg quantization allows the physical and PV photons to be handled in the same way, which 
facilitates the regularization and the preservation of symmetries. In Feynman gauge, this 
quantization is equivalent to the Gupta-Bleuler quantization used previously |3|, |5|, |8[. 

As a first application of the formalism, we have studied the dressed electron in a Fock 
space truncated to include at most one photon and no positrons. In particular, we have 
investigated the gauge invariance of the mass shift and the anomalous magnetic moment. 
In both cases, the residual gauge dependence can be ascribed to errors induced by the Fock- 
space truncation. The dependence of the anomalous moment on the gauge parameter C, is 
illustrated in Fig. SJin the case where the PV fermion mass mi is infinite. If the PV- fermion 
mass is kept finite, there is also gauge dependence of order 1/mi, due to fermion- flavor- 
mixing currents. The strong dependence as C is to be expected, because in this limit 
the gauge-fixing term is removed from the Lagrangian and the theory becomes undefined. 
We have also found that two PV-photon flavors are required to maintain the chiral symmetry 
of the massless-electron limit; this extends the earlier Feynman-gauge result (sl to arbitrary 
gauges. 

The formalism can be applied to higher-order truncations and to other charge sectors. For 
high-order truncations, where calculations are done numerically jsj, the PV- fermion mass 
mi must usually be kept finite; however, with mi large enough, the gauge-dependent effects 
should be small. 
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Appendix A: An integral identity 



From |5|, we have already that J = M^/q. Thus, to show that J + A J = M'^{Io + A/q), 
we need only consider the gauge-dependent parts. If we write the common denominators in 
(13. 9 p as the difference of two terms, we obtain 



A J - M^A/n 



327r^ 



dydk'\ 



M 



2 + fc2 H 

0^ 



(Al) 



M2 



y 



M2 



1-2/ 



r,2 



where we have used 

.2 



that contain m| + k]_, we change the integration variable y to 



— C)/^z^ to simplify the leading factors. For the terms 



or 



X 



X 



y) 



m]y + lifil - y) + kl 



2 ' 



(A2) 



(A3) 



depending on whether /i; or fii appears in the denominator of the integrand. As discussed 
in jsj, these variables range between and 1, though in reverse order, and satisfy 



"^j + f^±^l^f + kl_ + kl ^ + k 



y 



y 



X 



X 



and the analogous expression for x. With use of this identity and differentiation of ( ]A2I) . we 
obtain 

dy ^^(j^2_ + _ fif + kl 

' y{l — yY X \ 1 — X X 



(A4) 



(A5) 



The analogous expression holds for dx. Substitution into f lAll) and replacement of x and x 
by y, to have a common integration variable for all terms, yields 



A J - M^A/o 



327r2 ^ 



M2 



dydk 



(A6) 



X 



+ 



1 



M2 



M2 



y 



M2 



i-j/ 



y 



M2 



m^+k^ 

1-y 



3/ 

M2 



M2 



1-J/ 



2/ 
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Here we have taken into account the reversed order of hmits for x and x by changing the 
sign of the terms in the square brackets. Clearly, the sum of terms in the curly brackets is 
zero, and, therefore, A J = M^AJq. 



Appendix B: Evaluation of the anomalous- moment formula 



On substitution of the two-body wave functions (13.51) and use of the Kronecker deltas in 
spin, the expression fl3.27p for the anomalous moment becomes 



dydk] 



yn{P 



M2 



1-y 



(Bl) 



y 



X < e\>{k)V'jP -k,P){—+t, ^ ^ 



dk^ ' "dk^J j^p "^l+fej 



+e'CXmi'l{P-k,P) 



d . . d\ e\^\k)Vi:„_{P-k,P) 



dk^ dk"^ 



M2 



1-1/ y 



with the vertex functions specified in fl2.2ip and (12. 22 p . The terms generated by differen- 
tiation of the denominators cancel. Simplification of the remaining terms, summed over 
polarizations A, yields 



ji i'i" 



(B2) 



mi/(l -y) - rrij 



1/(1 -y) 



i-y 



+ 



iTLiii — rrij 2 + {1 — ()/( 



1-?/ 



+ 



rriiit — rrij f mjniiiy jij 



fj,f{l-y) \ l-y yC 

( 



X 



v 



i-y y 



^2 _ _ 



1-?/ 



The k^ integrals have double poles, if /iq < M — mo or /xq < ^ — ""^o! following the earlier 
convention in [sj, we define the integrals by 



f{x)dx 



lim — 



V 



f{x)dx 



-V 



f{x)dx 

{x — a)2 v-^o 2ri [_ J x — a — r] J x — a + r] 
With or without the pole, we find 

dkl y\i-yr 



(B3) 



J^2 _ '"l+fej _ /^f+fcj 



(B4) 



16 



These leave the anomalous moment in the form 



j7 j'i 

with Aoe a gauge-dependent part given by 



^a, = -^—^^MY,Y.^-^y'''^'^''ii4^A^^''-^j) I dyy{l-y) 



j7 



\m]y + - y) - M^y{l - y) ^^^^ 

mjUij^y^ + fi^{l-y)/C 1 

[mh + fij{l-y)- M^yil - y)][m']y + iq{l-y)- M^y{l - y)] j ' 



Except for the normalization factors zf, this expression is IR safe; the photon mass 
can be set to zero. We first take the limit mi — )■ oo, which, with use of z^' = z^ = Zt, implies 
J2ii~^y^i ~^ ^0 and ^•(— l)*mjZj — )■ Mzq. Next, the second PV-photon flavor is removed 
by the limit /X2 — ^ oo, in which ^2—^0 and — )■ 1, to obtain 

a -> ^Mz^ \-(-lY f dy yi^-ymy-2iM-mo)] 



and 



-M{M-mo)z^oJ2(-^y [ 

1=0 



C 2vr ^ ^^^"^^ ' J '\mly + ^^^{l-y)-MM^-y) 



yil-y)[moMy^ + fif{l-y)/C] 



[mly + fif{l-y)- M^y{l - y)][mly + - y) - M^y{l - y)] 



(B8) 



Appendix C: Evaluation of the normalization condition 

On substitution of the two-body wave functions (13.51) . the normalization condition (^Z 
becomes 



xe£^ (fc)e£) (fc) [l-^;;(P - fc, P)Vp^{P -k,P) + - k, P)U^,,^{P - k, P)_ . 
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To simplify the expression, we first add and subtract the A = term with the denominator 
replaced by the denominator of the A 7^ terms. We then have 



ji i'i" 



(C2) 



X 



+ 



M2 - 



1-3/ y 



( 



1 



^2 _ "'j+^i _ A?+fci 



On use of the vertex functions in fl2.2ip and fl2.22p and of = z~ = Zi, we obtain, 



1-^0 . ^g^2 



(C3) 



M2 



rri: — 2mj(mii + mj//)(l — y) + mj/mj"(l — y) + kj^ 
2 — ^ 



1 — ( f ruj {rriii + 



c 



1 



1 



1 



X 



1-3/ 1/ 



M2 - 



i-i/ 



^ /i^ ^ (ml + miimi't)k]_ 



{l-yY y3(^i_yy y3 y{l-y) 



y{^-y) y^{^-y). 



The kj_ integrals are defined by (1B3P when a pole is present. The integrals needed are given 
by », 



k I dk 



1-y 



and 



/l I 



1-j/ y 



V(l - ln[|m> + /.f (1 - - M^yll - (C4) 



2/2(1 - (m> + /if (1 - y) - M2y(l - y)) (C5) 



X 



(2 ln[|m|2/ + //^(l - - m2|/(1 - y)\] + l) 
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where we have dropped infinite terms that cancel in the final expressions, due to either 
'}2ii~^y^i = or the difference between the two denominators, (M^ — "^(^y'^ — ^' )^ and 
(M^ — '"j^^^ — ^i-ii)2^ Substitution of these integrals yields 



ji i'i" 



m]y + ii^il -y)- IVPy{l - y) 



2y{m'^j — 2mj{mii + mj//)(l — y) + mj'mj//(l 



(C6) 



^ ' m^imi, + mi.)y{l - y) + ^^'^ ^' 



C 



C 



y 



-2 ( y - ^1-^ ) ln[|m^y + /.^(i _ ^) _ M2y(l - 



C 



y 



+- 



1 



1 



m]y + - - Af2y(l - y) ni]y + ^2(1 - y) - M2|/(l - ?/) 



X 



-(m^y + _ ^) _ ^2y^^ _ ^)(2 ln[|m2y + (1 - - M^l - 



+- 



?/(m^ + miirriiii) + 



2/i?(l-l/) 



In 



m]y + fifil - y) - M^y{l - y) 



^]y + - - M^yil - y) 



Collecting terms with logarithms, we find that multipliers containing 1/y cancel. The re- 
maining terms containing 1/y are not singular in the full expression, which can be arranged 
explicitly by adding ^^^-f- to the curly brackets of (lC6p . This additional piece makes no 
contribution to the sum over /, because ^;(— l)'^z = 0. The resulting expression for the 
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normalization condition is 



l = zl- zl + 



167r2 



ji i'i" 
1 



j+i+i'+i'' 



iiZi>Zi» I dy 



^"jV + /^f (1 - 2/) - M2|/(l - y) 



(C7) 
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1-C 
C 
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1 
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x(l-y) [m^2-M2(l-y)] 



+— [y(m| + mi/mi//) + 2[mj2 - M^{1- y)]] In 
In the mi oo limit, this becomes 



m]y + il^^{l-y)-M^y{l-y) 



m]y + - y) - M2y(l - y) 
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47r 



dy 
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1 
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mly + - y) - M2y(l - y) ( mgy + ;i2^1 - y) - M2y(l - y) 
X (1-y) [mg-M2(l-y)] 
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